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Abstract 

We revisit the notion of approximate error correction of finite dimension codes. We derive 
a computationally simple lower bound on the worst case entanglement fidelity of a quantum 
code, when the truncated recovery map of Leung et. al. is rescaled. As an application, we 
apply our bound to construct a family of multi-error correcting amplitude damping codes that are 
permutation-invariant. This demonstrates an explicit example where the specific structure of the 
noisy channel allows code design out of the stabilizer formalism via purely algebraic means. 



PACS numbers: 03.67.Pp, 03.67.Lx, 03.67.Hk, 03.67.-a 



I. INTRODUCTION 



Quantum information, when left unprotected, often decoheres because of its inevitable 
interaction with the environment. The field of quantum error correction arose from the need 
to combat decoherence in quantum systems, and treats the decoherence as a noisy quantum 
channel. An important problem in quantum error correction is that of determining the utility 
of a given code with respect to the noisy quantum channel. The quantum error correction 
conditions of Knill and Lafiamme [1] are equations from which one can determine whether 
a quantum code is entirely robust against a given set of Kraus effects of the noisy channel. 
The Knill-Lafiamme conditions lie at the foundations of Gottesman's stabilizer formalism 
[12] from which quantum error correction codes are designed and studied. 

In this paper, we revisit the approximate error correction of finite dimension codes via a 
perturbation of the Knill-Lafiamme conditions. We derive a computationally simple lower 
bound on the worst case entanglement fidelity of a quantum code, when the truncated 
recovery map of Leung et. al. [2] is rescaled to guarantee its validity as a quantum operation. 
Our lower bound arises from repeated application of the Gershgorin circle theorem on the 
relevant matrices. 

The simplicity of our bound comes at a price - we do not have the near-optimal guarantees 
that the methods of Barnum-Knill [3] and Tyson-Beny-Oreshkov [1H6] yield. However in 
this trade-off, we are able to construct a family of multi-error correcting amplitude damping 
qubit codes that are permutation-invariant. We thereby demonstrate an example where the 
specific structure of the noisy channel allows code design out of the stabilizer formalism via 
purely algebraic means, as opposed to optimization techniques [714T0] and other approaches 
[T3HT5] . Our qubit permutation- invariant codes also extend the existing theory of qubit 
permutation- invariant codes pTl [T71 - [T9] : while no qubit permutation-invariant code corrects 
arbitrary single qubit errors, there exist qubit permutation-invariant codes that correct 
multiple amplitude damping errors. 

A. Organization 

In Section |Tl| we introduce notation and concepts needed for this paper, including 
quantum channels, quantum codes and the entanglement fidelity of a code. In Section 
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Ill we address the perturbed Knill-Laflamme conditions, revisit the Leung et. ah recovery 



map and determine using Lemma in.3 when the Leung et. aL recovery can be rescaled to 



a quantum operation. In Section III C , we prove our algebraic lower bound on the worst 



case entanglement fidelity (Theorem III.4). Finally in Section IV, we apply our lower bound 
to construct a family of qubit permutation-invariant codes that correct multiple amplitude 



damping errors (Theorem IV. 3). 



II. PRELIMINARIES 

A. Quantum channels 

For a complex separable Hilbert space let ^B(H) be the set of bounded linear operators 
mapping "H to "H. Define the set of quantum states on Hilbert space "H to be 25 ("H) where 
D{'H) is the set of all positive semi-definite and trace one operators in ^{H). For any 
subspace C G H with orthonormal basis Be, define 2)(C) to be the set of all elements in 
DiTi) that are invariant under conjugation by the projector 11 = J2\fs)eBc \ f^)(f^\- 

A quantum operation $ : ^{H) — ?■ ^{H') is a linear map that is completely positive 
and trace non-increasing. A quantum channel $ is a trace preserving quantum operation. 
In this paper, the Hilbert spaces H and H' are always isomorphic. A quantum operation 
$ : 03(7/) — i- 03(7^) can always be expressed in the Kraus representation \25\ : 

where p G ^(Ti), .^<i> C ^(Ti) is a set of Kraus operators of quantum operation $, and 1^ 
is the identity operator on complex Hilbert space H. 

B. Quantum codes and Entanglement Fidelity 

Define the minimum eigenvalue of a Hermitian matrix H restricted to subspace C to be 

Amin,c(H) := min (/3|H|/3). (1) 

|/3>ec 

The entanglement fidelity of a state p with respect to the quantum channel 7^ o ^ is 

Fe(p,7^o^)= ^ |Tr(Bp)|^ (2) 
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where the set of Kraus operators of 7^ is ^-jioA 126]. The entanglement fidehty of a state p 
with respect to the quantum channel TZo A quantifies how well the entanglement consistent 
with state p is preserved when the noisy channel is A and the recovery map is TZ. 

III. RESCALED TRUNCATED RECOVERY 

In this section, we analyze the performance of the rescaled truncated recovery map. 



Firstly in Section 111 A we study the diagonalization of the Knill-Lafiamme conditions. We 



next analyze the rescaling of the truncated recovery of Leung et. al. [2] into a quantum 



operation. In Section 111 C we prove Theorem 111.4, the main result of this paper. 



A. The perturbed Knill-Laflamme conditions 

In this subsection, we apply the canonical procedure [26j to diagonalize the error 
correction perturbed Knill-Lafiamme conditions. We also introduce notation that is used 



both in the statement and the proof of Theorem 111.4 We start by introducing the following 
notation. 

Nl. Let the M- dimension code C he a subspace of the Hilbert space H with projector U and 
orthonormal basis Bq- Let A : ^{Ti) — t- ^{T-L) be a quantum channel with truncated Kraus 
set Q. 

Using the notation of ]N[i| define an orthonormal basis {|E) : E G 17} C C'^' labeling the 
effects in Q. For all A, B G i7 and G Bq, define 

e(A,B,|a),|/3)) := {a\A^B\P) - gj^,BS\a),W). (3) 

to quantify the perturbation to the Knill-Lafiamme condition, where 

^A,B ■■=^J2 (/5|A^B|/3). (4) 



Define the Hermitian matrix 



M 

l/3>6Bc 



G = Yl ^?a,b|A)(B|. (5) 
A,Ben 
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The hermiticity of G implies the existence of a unitary matrix V and diagonal matrix D 
such that 



V= 5^ t;E,F|E)(F| (6) 
E,Fen 

D := VGV^ = ^rfE|E)(E| (7) 



which implies that 



y ] Ve,FV^'^F'9F,F' — C?E^E,E'- (8) 



F,F'en 

For all A G il, define the transformed Kraus operators 



A := 5^ t;A,FF. (9) 



F^n 

Substituting ([s]) into ^ gives 

(a|AtB|/3) = dA6A,B6ia),m + e(A,B, (10) 

where 

?(A,B,|a),|/3)) := (i;l,p,t;B,F)e(F', F, |/3)) . (11) 

F,F'6f^ 



Equation (10) gives the 'diagonalized' form of the perturbed Knill-Laflamme conditions. 



The transformed error is quantified by (11). Let £\a),\i3) '■= maxA.Bef^ e(^A,B, |a) 



Then the Cauchy-Schwarz inequality and normalization of the rows of V implies that 

|e(A,B, |/3))| < ^ \v*A^F''VB,F\e\a),\i3) < Me\a),\f3)- (12) 
F,F'eQ 

B. The Leung et. al. truncated recovery 

The truncated recovery of Leung et. al. j2j gives an algebraically simple lower bound on 
the entanglement fidelity of a code with respect to a noisy channel, under certain assumptions 
on the noisy channel and the quantum code. To understand this truncated recovery, we need 
the following notation. 
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N2. We use the notation 0/ Njlj For all A G define Ha '■= UaHUa^ where Ua is 
the unitary in the polar decomposition AH = UAV^nAAAn. Define the completely positive 
but not necessarily trace preserving linear operator 7ln,c '■ ^{'H) — )■ ^(T-L) where for all 
12 G m{H), 



:= Yl Ra/^RL Ra := UJ^Ha- 



(13) 



Aef^ 



Without loss of generality, we pick such that y/UA^AU is never the zero operator. 



We call the completely positive map TZn^ a truncated recovery map, because is a 
truncated Kraus set of the quantum channel A. The truncated recovery map TZq^ is also a 
quantum operation when the projectors Ha in ]N[2] are orthogonal, that is 



nAnB = nA5A,B VA,BGfi. 



(14) 



Lemma III.l (Leung et. al. [2J). In the notation 0/ ]N[T] and ]N[2| for any p G 2)(C), the 
bound J^Aen |Tr(RAAp)|^ > Z^Aen '^min,c(Al'A) holds. 

Proof. For A G fi, define ^a '■= Amin,c(A^A) and the positive semidefinite residue operator 

tta := VnAtAn - y^U. 

Substituting tta into the polar decomposition of AH gives 



An = UA^nAtAn 

= UA(VnAtAn - + VT^in) 

= UA(vrA + ^/JlAJ^) 

= UA(7rA + y/J^tn )n. 



(15) 



The spectral decomposition of p in the basis Be and equation (15) imply that 
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j2\n^AAp)f>j2 



E 

Aen 



Aen 



J2 P|/3>(/3|RaA| 

l/3>GBc 



J2 m(/3|(nul)(UA(v^i« + 7rA)n)|/3) 

l/3> eBc 



(16) 
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Using n|/3) = 1/3) and U^Ua = ^ becomes 

E( E m(/3i(v^%+^A)i/3))'. 

Aen |/3>eBc 

Moreover tta is positive semi-definite, hence the above expression is at least 

2 



Aen \i3)eBc Aef7 



□ 



Leung et. al. proved that when the orthogonahty condition (14) holds, the truncated 
recovery map TZn,c is also quantum operation, and thus Lemma III.l gives a lower bound 
on the worst case entanglement fidelity. We detail this in Lemma III.2 

Lemma III. 2 (Leung et.al. [2]). We use the notation o/ Njl] and ]N[2| Suppose that (I4) 
holds. Then the truncated recovery map TZn,c is a quantum operation and 



min Fe(p, 7^n,c o A) > V] A^in,c(AtA). (17) 



Proof. Observe that 



J2 Ri^A = Yl nAUAU^HA = J2 ^A- (18) 

Aen Aen Aen 

Orthogonality of the projectors IIa implies that 

E "a ) ( E nA' ) = Yl nA<5A,A' = E "a. (19) 

^Aen / \A'en / A,A'en Aen 
Hence XIagh R-aR-a is also a projector. The map 7ln,c is also completely positive, that is 
'^Aen I^ApR-A — 0) because U^pUa > which implies that RapRa ~ nU^pUAH > 0. 
Obmitting non-negative terms in the sum pertaining to the entanglement fidelity, we get 

J'e{p,nn,c o ^) > ^ I Tr(RAAp)|2. (20) 

Aen 



Applying Lemma III.l gives the result. □ 



Rescaled maps have been used in the study of near-optimal quantum recovery operations, 
including the Barnum-Knill recovery map f3] and the Tyson-Beny-Oreshkov quadratic 
recovery map [H [5] . In the notation of INjl] and ]N[2| the completely positive map TZn^ 
might increase trace and hence not be a quantum operation. Fortunately a bounded TZn,c 
can be rescaled to the quantum operation TZnfi^rj '■= (1 + v)^^ T^nfi- 



Lemma III. 3. Using the notation 0/ INljlj and ]Nj2| let rj > maxA^Bef 
Then 7ln,c,rj is a quantum operation. 



nu^Uen 



Proof. It suffices to show that | J^Aen R-a^^a < 1 + First observe that 

J2 RLRa = Yl nAUAU^HA = Yl ^A- 

Aen Aen Aen 

The projectors 11 a may not be orthogonal, so 



\Aen 



Aen 



Ben 



A,Bef^ 



:^nA+ Y UAnuiUgnu 

Aen Af^Ben 



and hence 

Y<^ 

'ash 

Let e = II^R^Ra 



J^RLRa llE^A^i 



+ max 
2 A^Ben 



nu^Uen 



Aen 



Aen Aen 
1 > 0. Then substituting e and f] into the above inequahty gives 



l + 2e + e^ <l + e + r], 



which imphes that e < rj. Hence 
T^n^Cri is a quantum operation. 



T^Aen R-A^-^ 



< l + rj, and the completely positive map 

□ 



C. Algebraic lower bounds on the worst case entanglement fidelity 

Here, we prove algebraic lower bounds on the worst case entanglement fidelity of a code 
using a rescaled truncated recovery map, given partial knowledge of the noisy quantum 
channel. Our main technical result is the following: 

Theorem IIL4. Let the M-dimension code C with an orthonormal basis Be be a subspace 
of the Hilbert space Ti. Let A : 03(7/) 03(7/) be a quantum channel with truncated Kraus 

set n. Define G = Y.A.Ben9A,B where gA,B ■= Xl|/3>eBc ^^l"^^^'^)' suppose that 
Amin(G) > 0. Suppose that for all Kraus effects A,B in the truncated Kraus set Q and for 
all distinct orthonormal basis vectors \a), \f3) G Be, 



\9a 



,B - (a|AtB|a)| < e, I (a|A"^B|/3) I < e. 



(21) 
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Then the minimum entanglement fidelity of our code C with respect to the noisy channel 

13, N -1 



A zs at least (Tr G - M\n\^e) 1 



M\n\-^e 

AmiTi(G) 



The proof of Theorem III.4 follows from the direct application of Lemma IIL3 which gives 
us the important properties of our rescaled truncated recovery, and the repeated application 
of the Gershgorin circle theorem. 



Proof of Theorem III. 4 From ([s]), we have 



|a>eBc 



which implies that 



^ e(^A,B,\a),\a))=0. 



\a)el3c 



Hence Amax(nA"l'An) > ^ ^ (a;|A^A|a) > c^a which is at least Amin(G)- 



\a)£Bc 



Applying the Gershgorin circle theorem on the matrix nA'''Bn with the error estimate 



(12), we get Amax(nAtBn) < M\Q\e for distinct A, B G fi. 



For distinct A,B G fi, let All and BII have polar decompositions AH = U^vnAtAII 
and BH = Us VnBtBn respectively. Then 



UA^BU = V HAt Anut Ug V HBtBH = V HAt An(nut UgH) V HBtBH. 
Hence by the sub-multiplicative property for norms, ||nU~^Ugn||2 is at most 



HAtAH 



Hence by Lemma 



HA+BH 



VnBt 



BH 



< 



iHAtBHl 



< 



M\n\e 



min Amax(nF'fFn) A„iin(G) 

Fe{A,B} 



in.3 



the map TZn,c,ri is a quantum operation whenever 77 > 
the Gershgorin circle theorem, Amin,c(A^A) is at least c/a — M|f2|e, and hence 



'^min(G) 



By 



Y Amin,c(A"^A) > Y^^^ ~ M\n\e) = TrD - M\n\'^e = TrG - M|npe. 



Aef^ 



Use of Lemma III.l then gives the result. 



□ 



Under stronger assumptions on the set of conditions the truncated set of our noisy channel 
must satisfy, we can obtain the following corollary. 
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Corollary III. 5. Let t be a positive integer, and 7 > 6e an error parameter, and assume 



that the requirements of Theorem II 1. 4 hold with e = 'j'^^'^^. Suppose that every A E Q is of 
order 0(7''^) for some non-negative integer Sa < t, and all Kraus effects not in Vt are of 
order (9(7*"*^^). Then the minimum entanglement fidelity of the code C with respect to the 
noisy channel A is at least 1 — 0(7*+-'^). 

Proof. Note that Tr G = 1 - 0(7*+^), and X^^{G) = 0(7*) ^ 0. Substitution of these 



parameters into Theorem III.4 gives the result. □ 



IV. APPLICATION: PERMUTATION-INVARIANT AMPLITUDE DAMPING 
CODES 



In this section, we apply Theorem III. 4 to prove the existence of family of multiple 



amplitude error correcting permutation-invariant codes. 

First we define permutation-invariant codes encoding a single qubit. Define Pm.a to 

be the set of all length m binary vectors of weight a, and the corresponding permutation- 

/ \ -1/2 

invariant qubit states to be \Pm,a) '■= XlxePm a l-^M ) " "^^^ permutation- invariant codes 
we consider have the logical codewords IJl) := Yl^=o ^/■^j,a\P'm,a) where j G {0, 1} and IJl) 
are linearly independent. Moreover the non- negative weights Xj^a are normalized so that 

Now we introduce notation relevant to the amplitude damping channel. The amplitude 
damping channel has Kraus operators Aq = |0)(0| + ^/T^^\l){l\ and Ai = y7|0)(l|, 
where the non-negative parameter 7 < 1 quantifies the amount of amplitude damping of 
Ay. We denote the Kraus operators of by A^ := A^-^ ... ® A^^, with k = {ki, km) 
is a binary vector. For positive integers t, we say that a code C is a t-amplitude damping 
code (or a t-AD code) if there exists a quantum operation IZ such that 

max^(l-F,(p,7^o^^™)) =0(7*). (22) 

The recovery operation that we use in this section is the rescaled recovery map of Leung et. 
al. 
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Let k and k' be length m binary vectors with weights k and k' respectively. Then 

m 



6=0 



b-k 



\ \b + k'-k 



m ~ T 



b-k 



(23) 



where r is cardinality of the union of supports of k and k'. The expression (23) allows us 



to deduce algebraically sufficient conditions on A, a for our code to correct t-AD errors. 



Lemma IV. 1. Let m and t be positive integers, with m > t. Let Xj^h = for all integers 
6 G (m — t,m] and j G {0, 1}. Further suppose that for all non-negative integers c < t and 
£ < 2t, we have 



m — b\ fb 



m — b\ fb 



fe=0 ^ / \ / b=0 

Let k and k' be binary vectors of equal weight k where k <t. Then 



(0L|4^k'|0L) = (iLlAUk'IU) + 0(7'*+'). 



Proof. For non-negative integers m, fc, 6, c such that b < m — c and c < k, we have 

m — k — c\ /b\ 

b — k I \k) fra — b\ f k\ . 

dr. 



m 
k — c 



c / \ c 



Using the identity ( 24 ) with ( 23 ) , we get 



{jL\AiA^'\jL) = J2hbi'a-i) 



6=0 



m \ \ c J \c ' 
k-c) 

b 



6=0 /3=0 ^ 

k 



C! 



k ~ c 



c 

c / \k — c 



k — 



m 



6=0 

-1 m 

7 7 ,^i,6 

6=0 



/3=0 ^ 
b-k 



c / ^ \k + BJ \ k 

13=0 ^ ' \ 



(24) 



(25) 
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Hence the coefficient of 7^^+*^ in the expression above is 

:)<^o^(;!j'(gvC"-)G:,))cr)<-)' 

By the assumption of the lemma, the desired result follows. □ 

We also need the following combinatorial lemma: 

Lemma IV. 2. Let t be a positive integer. Then for any non-negative a <t + 1 we have 

t+i 



Proof. Let g{x) := (1 - so that g^'^\x) = (t + l)(a)(l - and ^^"^1) = 0. 

Substituting x = 1 into the expansion 



i=0 

we get the binomial identity 

t+i 



j=0 



We prove our lemma by inducting on a. The binomial identity (26) implies that our lemma 



is true for the base cases a = 0, 1. The identity (26) also implies that 



j=0 ^ ' /3=0 \i=0 ^ 



Expanding the falling factorial ^(cl'+l) into a sum of monomials in % when a' < t, we get 

't^\\ 

for some choice of constants hp G R. Note that the bracketed term in the equation above 
is zero by the hypothesis that our lemma is true for a = a', where 1 < a' < t. Hence the 
lemma also holds for a = a' + 1. □ 

The existence of permutation-invariant amplitude damping codes is given by the following 
theorem: 
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Theorem IV. 3. Let t be any positive integer and m = 9t^ + 4t. For all j G {0, 1} and 

integers h <m divisible by 3t let 

1 + 



A 



■3,b 



\b/mj- 



,-3t ■ 



(27) 



and all other values of Xj^b be zero. Then the span of \0l) = YlT=o ^/^o,b\Pm,b) and = 
SfcLo \/'^i,b\Pm,b) is a t-AD permutation-invariant code. 

Proof. By definition, Aj^b satisfy tlie normalization condition Xlblo = 1, so Xj^b define 
valid logical codewords Ijl). Moreover, |0l) and have distinct support, and are hence 
linearly independent spanning a two-dimension codespace. 

Now define the truncated Kraus set of amplitude damping effects Vl := {A^ : k G 



Z^, ||k||i < t}, so that Vt satisfies the order constraints of Corollary III. 5 Define the matrix 
Gj to be ^E,FenOi|E"^F|jL)|E)(F|, and G := ^fl±^. To use Corollary |lIL5| we first have 
to prove that the matrix G is positive definite. 

Let A = Xlpen'^A.FF and V := XIe Fen '"e,f|E)(F| be as defined in (9) and 
(j6| respectively, and correspondingly define the vector |^a) 
serve that AminfG) = Ar 



EF6f1^A,F|F). Ob- 

^(VGVt) = minAen(A|VGVt|A) = minAen(*A|G|^A) = 



minAen^(^A|(Go + Gi)|^a) = miuAen ^ E]=oO'il = miuAen ^ Ej=o II A|jl)||2. 
Now the Kraus elements in VL annihilate at most t excitations, but the logical states Ijl) are 
permutation-invariant with support containing at least St excitations. Hence ||A|jL)||2 > 
which implies that Amin(G) > 0. 

Our choice of Xj^ implies that (ji^lAj^Ak'ljL) = and (OL|Aj^Ak/|li) = when k, k' G 



such that ||k||i = ||k'||i < t. Now we set t to 3t in Lemma IV. 2, and note that the coefficient 



of Xj^b in the assumption of Lemma IV. 1 is a polynomial of order no more than 3t in the 



variable h. Then applying Lemma IV. 2 to Lemma IV. 1 shows that the conditions of Theorem 
III. 4| are satisfied with e = 0(7^*). □ 
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